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ABSTRACT 

A luminosity-temperature relation for clusters of galaxies is derived. The 
two models used, take into account the angular momentum acquisition by the 
proto-structures during their expansion and collapse. The first one is a modifi- 
cation of the self-similar model (SSM) while the second one is a modification of 
the Punctuated Equilibria Model (Cavaliere et al. 1999). In both models the 
mass-temperature relation (M-T) used is based on the calculations of Del Popolo 
(2002b). 

We show that the above models lead, in X-rays, to a luminosity-temperature 
relation that scales as L oc T^, at scale of groups, flattening to L oc for rich 
clusters and converging to L oc at higher temperatures. However a funda- 
mental result of our paper is that the non-similarity in the L-T relation, can be 
explained by a simple model that takes into account the amount of the angular 
momentum of a proto-structure. This result is in disagreement with the widely 
accepted idea that the above non-similarity is due to non-gravitating processes 
as those of heating/cooling. 

Subject headings: cosmology: theory - large scale structure of universe - galaxies: 
formation 



1. Introduction 

Observations of clusters of galaxies (e.g, ROSAT, ASCA), performed in the past decade, 
have shown the existence of a tight correlation between the total gravitating mass of clus- 
ters, Mtot, their X-ray luminosity (Lx), and temperature (Tx) of the intra-cluster medium 
(ICM) (David et al. 1993; Markevitch 1998; Horner, Mushotzky & Sharf 1999 (hereafter 
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HMS)). The importance of these relations is due to the fact that cluster masses are difficult 
to measure directly, and when comparing cluster observations with models of structure for- 
mation a surrogate for cluster mass is used. Since M^ot compares with the ICM temperature 
measurements which can be obtained through X-ray spectroscopy, this explains the impor- 
tance of an M-T relation. On the one hand, the X-ray temperature measures the depth 
of the potential wells, and the bolometric luminosity, L oc n'^R^T^^'^, emitted as thermal 
bremsstrahlung by intra-cluster plasma measures the baryon number density, n within the 
volume i?x. Till some years ago, the cluster structure was considered to be scale-free, which 
means that the global properties of clusters, such as halo mass, luminosity-temperature, and 
X-ray luminosity would scale self-similarly (Kaiser 1986). In particular, the gas temperature 
would scale with cluster mass as T oc M^/'^ and the bolometric X-ray luminosity would scale 
with temperature as L oc T^, in the bremsstrahlung-dominated regime above ~ 2 keV. 

Studies following that of Kaiser (1986) showed that the observed luminosity-temperature 
relation is closer to L oc (e.g.. Edge & Stewart 1991), indicating that non-gravitational 
processes should influence the density structure of a cluster's core, where most of the lu- 
minosity is generated (Kaiser 1991; Evrard & Henry 1991; Navarro, Frenk & White 1995; 
Bryan & Norman 1998). One way to obtain a scaling law closer to the observational one is 
to have non-gravitational energy injected into intra-cluster medium (ICM) before or during 
cluster formation. This solution, called pre-heating, was originally invoked to solve two re- 
lated problems: a) to explain (Kaiser 1991; Evrard & Henry 1991) ^ the apparent negative 
evolution of the X-ray cluster luminosity function (Gioia et al. 1990; Henry et al. 1992) from 
the Einstein Medium Sensitivity Survey in a ilm = 1 Universe ; b) to explain (White 1991) 
^ why groups and low-mass clusters seem to have higher X-ray temperatures than expected 
based on member velocity dispersions. 

The mechanisms proposed to explain the slope change of the L-T relation 
can be divided into three main categories: 

(i) models that include a pre-heating of the gas within a cluster. Ponman et al. 

(1999) showed that the entropy of the ICM in the centre of low-temperature clusters is greater 
than the value expected from gravitational collapse. It has been shown that models that 
include an additional gas entropy can successfully reproduce many observational properties 



''indeed, numerical simulations that include gas dynamics but exclude non-gravitational processes such as 
radiative cooling and supernova heating produce clusters that obey these scaling laws (e.g., Evrard, Metzler, 
& Navarro 1996; Bryan & Norman 1998; Thomas et. al 2001b). 

^Kaiser's self similar model predicts L oc T^-^. Evrard & Henry (1991) obtained the relation L oc T^^/^ 

^In this case pre-heating was in form of supernovae- driven galactic winds 



-3- 



(Bower et al. 1997; Cavaliere et al. 1997, 1999; Tozzi & Norman 2001; Borgani et al. 2001; 

Voit & Brian 2001). 

(ii) models that implement feed-back processes that alter the gas characteristics 
during the evolution of the cluster. In principle, there are many different physical 
processes that could break the self-similar scaling, including heating from SN or from AGN, 
or the removal of low-entropy gas via radiative cooling with subsequent supernova heating 
(Voit & Bryan 2001). Another example is that of Muanwong et al. (2001), who simulated 
galaxy cluster formation including radiative cooling with cool gas dropout and was able to 
reproduce the L ocT^ dependence, without adding any entropy to the gas. Moreover other 
possibilities, such as magnetic pressure or cosmic-ray pressure have not been ruled out. Allen 
& Fabian (1998) have examined the effects of cooling flows for a sample of the most X-ray 
luminous (I/boI > 10'^^ erg/s) finding a flattening from L oc T"^ to L oc in agreement 
with models that include the effects of shocks and pre-heating on the X-ray gas (Cavaliere 
et al. 1997, 1999). Cavaliere et al. (1997, 1998), have constructed a model in which the 
observed L-M relation on both cluster and group scales can be reproduced by varying the 
gas density at the virial radius, according to the accretion-shock strength, as determined by 
the temperature difference between the infalling and virialised gases. Another possibility to 
explain the L-T relation are systematic variations in the baryonic fraction with cluster mass 
(David et al. 1993). To distinguish among these processes, observations of high-redshift 
groups and clusters will be crucial to measure the evolution of the observed 
scaling relations as function of redshift . 

(iii) hydro- dynamical models that do not include gas pre-heating, nor feed-back 
processes, which also reproduce the available observational data (e.g., Bryan 8z 
Norman 1998). Throughout this paper, we will analyse this last scenario. 

On the other hand, the mass-temperature relation, seemed like it ought to be more 
fundamental and less sensitive to non-gravitational effects. Yet, observations collected over 
the last few years indicate that this relation also disagrees with both the scale-free predic- 
tions and simulations that exclude non-gravitational processes (Horner, Mushotzky, & Scharf 
1999; Nevalainen, Markevitch, & Forman 2000 (hereafter NMF); Finoguenov, Reiprich, & 
Bohringer 2001 (hereafter FRB); Xu, Jin, & Wu 2001). These results derive mostly from 
resolved X-ray and temperature profiles coupled with the assumption of hydrostatic equilib- 
rium, and they do seem consistent with gravitational lensing measurements (Allen, Schmidt, 
& Fabian 2001). Understanding the scaling properties of clusters is of broad importance be- 
cause these scaling laws are integral to determination of cosmological parameters. Thus, any 
inaccuracies in the mass-temperature relation propagate into uncertainties in cosmological 
parameters derived from clusters (e.g., Voit 2000, hereafter V2000). 

In Del Popolo (2002b), we derived the mass-temperature relation and its time evolution 
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for clusters of galaxies in different cosmologies. We use two different models: the first one is 
a modification and improvement of a model by Del Popolo & Gambera (1999), based upon 
a modification of the top-hat model in order to account for angular momentum acquisition 
by proto-structures and for an external pressure term in the virial theorem. The second one 
is an improvement of a model proposed by V2000, again to account for the angular 
momentum acquired by proto-structures during their formation. Both models 
showed that the M-T relation is not self-similar. A break is present in the quoted relation at 
T ~ 3keV and, at the lower mass end, the power law index of the M-T relation is larger than 
q; = 3/2 even in flat universes. The slope of the power- law index depends on the considered 
cosmology. The two models also agree in predicting a more modest time evolution of the 
quoted relation in comparison with the results of previous models, which also depends on 
the cosmology. 

This is in agreement with studies showing that the self-similarity in the M-T relation seems 
to break at some keV (NMF; Xu, Jin & Wu 2001). By means of ASCA data, using a small 
sample of 9 clusters (6 at 4 keV and 3 at ~ 1 keV), NMF has shown that Mtot oc T^'^^^^-^^ for 
the whole sample, and Mtot oc T^"^ excluding the low-temperature clusters. Xu, Jin & Wu 
(2001) has found Mtot oc T^-60±o-04 ^^g^^g p model), and Mtot oc T^-8i±oi4 by means of 
the Navarro, Prenk & White (1995) profile. FRB have investigated the T-M relation in the 
low-mass end finding that M oc T~^, and M oc T~^/^ at the high mass end. This behaviour 
has been attributed to the effect of the formation redshift (FRB) (but see Mathiesen 2001 
for a different point of view), or to cooling processes (Muanwong et al. 2001) and heating 
(Bialek, Evrard & Mohr 2000). Afshordi & Gen (2001) (hereafter AG) have shown that non- 
sphericity introduces an asymmetric, mass dependent, scatter for the M-T relation altering 
its slope at the low mass end (T ~ 3 keV) . 

The L-T and M-T relations are somehow related: as shown by Shimizu et al. (2003), it is 
possible to make a reliable prediction for the L-T relation once the M-T relation is specified. 
In turn, one can obtain the M-T relation that reproduces the observed L-T relation without 
assuming an ad hoc model for the thermal evolution of intra-cluster gas. The two relations 
(M-T and L-T) are strictly connected, as shown by Shimizu et. al (2003). 

This conclusion in turn indicates that the L-T relation provides a good diagnosis of the 
underlying M-T relation, which is as yet poorly determined observationally. 

Apart from the physical mechanism of the additional thermal processes, there are three 
effects that might modify the mass dependence of X-ray luminosity and steepen the resulting 
L-T relation. First, the gas density profile might be significantly fiatter for less massive 

1/3 

systems. Second, the mass dependence of the hot gas mass fraction is strong as /gas oc M^/, . 

2/5 

Finally, the mass-temperature relation is Tgas oc M^^ . In practice, a realistic model should 
be a combination of those three effects to some extent. 
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In this paper we derive a luminosity-temperature relation for clusters of galaxies that 
takes into account the amount of the angular momentum of proto-structures. We use two 
different models: the first (that we call MSSM) is a modification of the self-similar model 
(SSM) while the second one is a modification of the Punctuated Equilibria Model (hereafter 
MPEM) (Cavaliere et al. 1999). We show that the presence of the angular momentum during 
the gravitational collapse leads to non-self similar L-T relation. The two models used are 
described in Sect. 2. The results are presented and discussed in Sect. 3 and the conclusions 
are summarised in Sect. 4. 

2. Model 

2.1. Modified Self-Similar model for the L-T relation 

The L-T relation constitutes a fundamental hnk between the physics of the baryon 
component and the dynamical properties of the Dark Matter condensations. The simplest 
model describing that relation is the SSM model (Kaiser 1986), obtained assuming that the 
gas density or the baryon number density, n, is proportional to the average Dark Matter 
density, p, and that the virial radius, -Rvir is proportional to Rx (see introduction for a 
definition). In this way, one obtains, according to this last, L oc M^i^pT^^"^. In fact, T oc 
Mvir/i?vir, n oc p oc Mvir/i?^i,, i?vir OC Rx and L oc j^"'" p^T^''^r'^dr oc p'^T^/^R^^^, and 
recaUing that i?vir oc (Mvir/p)^/^, leads to L oc pM^i^-T^/^ or recalling that i?vir oc {T/pY^^, 
we get L oc p^/^T^. This last result is inconsistent with observed correlation close to L oc 

(Edge & Stewart 1991; Mushotzky 1994). Additionally, a further steepening at the 
temperature of galaxy groups is indicated for thermal emission not associated with single 
galaxies (Ponman et al. 1996). 

In the following, we derive a modified SSM, showing that slope of the L-T relation 
changes at different scales. 

Using Balogh et al (1999) notation, let us begin with a cluster with gas temperature 
T{r), density profile p{r)g, for which the bolometric X-ray luminosity from Bremsstrahlung 
scales as: 



(see Balogh et al. 1999), where Ci = 3.88 x lO^^s K^^/^ ^^^^-3^ ^ ^ o.59, 7?vir oc {M^i,/pf/^ is 
the virial radius where p{z) oc (1 + z)'^ is the Dark Matter (DM) density in the cluster, pro- 
portional to the average cosmic DM density pu at formation. The simplest model describing 
the L-T relation, that can be calculated by Eq. (1), is the SSM (Kaiser 1986), assuming that 




(1) 
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the pg (X p. 

We only consider haloes in which not all of the gas within ii^ir has had time to cool 
since the halo formed. 

We assume a singular, truncated isothermal sphere for the dark matter potential, p(r) = 
Pij(r/i?vir)^^, where pr is the density at the virial radius, i?vir, and is equal to a third of 
the mean density within R^ir, p(-Rvir)- This latter quantity is related to the critical density 
at redshift z by p(i?vir) = Ac{z)pc{z), and = 78Q(z) + 80 + 300Q(z)/(l + 15Q(^)) is 
a fit, accurate to better than 2 per cent, to the results of the spherical collapse model as 
presented in Eke et al. (1996). It will be convenient to define a redshift evolution term, 
Fi{zf = (1 + zf{l + Oo-2)Ae(^)/Ae(0), SO that 

PR = ^A,{0)p,{0)F,{zf. (2) 

For flo = 1, Fi{zY = (1 + 2;)^ and Ac = 178. In this model, we make the common assumption 
(e.g.. Eke et al. 1996) that the gas is distributed isothermally, with a temperature equal to 
the virial temperature of the halo. If the gas is dissipationless, its density profile will match 
that of the dark matter, i.e., 

Pg{r) = pg,R{r/R^ir)~^, (3) 

and Pg^B./pR = Qb/^o- To avoid the singularity at r = when integrating over the assumed 
isothermal profile, an arbitrary core radius of Vc = /c-Rvir is adopted with fc = 0.1, such that 
Pg{r < Tc) = Pg{rc)- The integral in Eq. (1) is dominated by the contribution from within a 
few core radii, and thus the scaling properties of this integral depend weakly on the assumed 
density profile. Furthermore, departures from the standard profile can be accommodated by 
redefining the core radius of the system. 



In order to obtain the luminosity-mass relation, we evaluate Eq. (1) and use the M-T 
relation, which takes into account the angular momentum of the proto-structure, obtained 
in Appendix A: 



kT ~ 8keV 



Mi 



1 I (kl\i I Ki{mi,x) 
mi ^ \ t } ^1 M8/3 
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3 _l_ Ko{mi,x) 



(4) 



where Ki{mi,x) ^ is given by: 

Ki{mi, x) — (mi — 1) FxLerchPhi{x, 1, 3mi/5 + 1) 



Ko{mi,x) indicates that Ki{mi,x) must be calculated assuming t = to 
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(mi - 1) FLerchPhi{x, 1, 3mi/5) (5) 

and 

mi = 5/(n + 3), to = ^ 3-, x = 1 + (i^V/^ which is connected to mass by M = 

i^o(i-f2o-nA)^ ' * 

MoX"^™'^/^ (V2000) and = where r^a is the turn-around radius and Xi is defined by the 
relation M — with ph is the background density Finally we get: 



3.31 X 10"^ ( 4^ I 



178 



mi ^ V t / /jwNS/s 



1-/C 
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which differently from Kaiser's (1986) prediction is not self-similar. It reduces to the self- 
similar form (L oc M^/^) if angular momentum acquisition is not taken into account, namely 
if £ ^ (or F ^ 0). 

The previous computation depends on the value of the angular momentum acquired by 
the DM haloes from tidal torques from surrounding matter. This enters the L-T relation 
through the quantities F and Ki (see also Appendix A). In the limit of vanishing angular 
momentum the L-T and the M-T relations reduce to the well-known self- similar forms. Then, 
it is important to add a discussion on the magnitude of the angular momentum calculated 
as in previous papers (e.g., Del Popolo & Gambera 1998; Del Popolo et al. 2001). 

The angular momentum is acquired by the cosmological torque acting on the proto- 
structures due to the tidal field of the environment. The amount of angular momentum 
as well as its distribution are related to the assumed power spectrum of density perturba- 
tions. We have to note here that the problem of the growth of angular momentum of proto- 
structures from the tidal torques of the surrounding matter has been studied extensively 
in the literature with both analytical and numerical (N-body) methods (e.g. Efstathiou & 
Jones 1979; Barnes & Efstathiou 1987; Voglis & Hiotelis 1989; Warren et al. 1992; Eisenstein 
& Loeb 1995; Kratsov et al. 1998). A main result of the above studies is that the values of 
the dimensionless spin parameter A = C\E\^^'^/GM^^'^, (Peebles 1971), follow a log-normal 
distribution with a small average value 0.05. In the above relation C is the total angular 
momentum of the proto-structure, E is its binding energy, M its mass and G the gravita- 
tional constant. The above numerical results are confirmed by analytical studies presented 
by other authors as those of Steinmetz & Bartelmann (1995) and Catelan & Theuns (1996). 
To be more precise, A depends on the galactic morphological type, being as high as A ~ 0.5 
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for spirals and SO galaxies, and A ~ 0.05 for ellipticals, although the dispersion around these 
values is large (Efstathiou & Jones 1979). In the case of structures of 10^^ — 10^^ its value 
is ~ 0.1 and ~ 0.01 for clusters ^. 

In this paper, we calculated angular momentum as Sect. 3 in Del Popolo et al. (2001), 
following Eisenstein Sz Loeb (1995). With the Bardeen et al. (1986) power spectrum 
smoothed on galactic scale for a z/ = 2 peak, the model gives a value of 2.5 x lO^^gcm^/s, in 
very good agreement with Catelan & Theuns (1996), in other words the amount of angular 
momentum used in our calculations is consistent with the values of A predicted by the tidal 
fields of the surrounding matter. Although this amount is in general small, our results show 
that it is efficient to lead to a non-similar L-T relation. The angular momentum of dark 
matter haloes has also other important consequences. For example small amounts 
of angular momentum are able to change the density profile of dark matter haloes from the 
isothermal law p(r) oc to a profile that it flattens significantly inwards, (e.g. Hiotelis 
2002). 

Moreover, several studies have shown that the influence and the role of shear on structure 
formation is of fundamental importance. Shear on a density perturbation can be produced 
by the intrinsic asphericity of the perturbation itself (internal shear) or it can be due to the 
interaction of the perturbation with the neighbouring ones (external shear). For example, 
according to the previrialization conjecture (Peebles & Groth 1976, Davis & Peebles 1977, 
Peebles 1990), initial asphericities and tidal interactions between neighbouring density fluc- 
tuations induce signiflcant non-radial motions which oppose the collapse. This means that 
virialized clumps form later, with respect to the predictions of the linear perturbation the- 
ory or the spherical collapse model, and that the initial density contrast, needed to obtain 
a given final density contrast, must be larger than that for an isolated spherical fluctuation. 
This kind of conclusion was supported by Barrow & Silk (1981), Szalay & Silk (1983), Vil- 
lumsen & Davis (1986), Bond & Myers (1993a,b) and Lokas et al. (1996). Arguments based 
on a numerical least-action method lead Peebles (1990) to the conclusion that irregularities 
in the mass distribution, together with external tides, induce non-radial motions that slow 
down the collapse. In a more recent paper. Audit et al. (1997) they conclude that spherical 
collapse is the fastest. This result is in agreement with Peebles (1990), and more recent 
papers, namely Del Popolo et al. (2001), Del Popolo (2002a). 



^The resulting typical circular velocities of structures is ~ 150km/s for galaxies similar to the Milky Way, 
~ 5 km/s for clusters and ~ 10 km/s for superclusters (see Catelan & Theuns 1996) 
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2.2. Improvements to the punctuated equilibria model 

In this subsection, we shall extend the punctuated equilibria model (PEM) by Cava- 
liere, Menci e Tozzi (1997,1998, 1999) (hereafter CMT97, 98, 99) to take account of angular 
momentum acquisition from the proto-structure. In their model (CMT98), the cluster evo- 
lution is described as a sequence of "punctuated equilibria" (PE), that is to say, a sequence 
of hierarchical merging episodes of the DM halos, associated in the ICP to shocks of various 
strengths (depending on the mass ratio of the merging clumps) , which provide the boundary 
conditions for the ICP to readjust to a new hydrostatic equilibrium. 

The X-ray bolometric luminosity of a cluster is given by Eq. (1), which in CMT98 
notation is: 

Loc / n'(r) TV2(r) . (8) 

Jo 

Here T{r) is temperature in the plasma and r2 is the cluster boundary, that we take to be 
close to the virial radius i?vir oc M^/^p^^/^, where p{z) oc (1 + zY is the DM density in the 
cluster, proportional to the average cosmic DM density Pu{,z) at formation. 

As shown in Appendix B, the L-T relation can be casted in the form: 

1/2 



L OC . 





"T2" 




[Tv\ 



\8/3 



mi ^ \ t I ^ 



\ 1 I /tn^2/3 Kq 

See Appendix B for a derivation of Eq. (9) and a definition of the terms involved. 



(9) 



Our final aim is to compute the average value of L and its dispersion, associated with 
a given cluster mass m. 

In order to reach this goal, we must sum over the shocks produced at a time t' < t 
in all possible progenitors m' (weighting with their number) by the accreted clumps Am 
(weighting with their merging rate); finally, we integrate over times t' from an effective lower 
limit t - At. 

The average L is then given by 

{L) = Q / dt' / dm' / dAm -^(m' ,t'm,t) — ^L; 10 

Jt-At Jo Jo dm'^ ' ' ' ^ dAmdt' ^ ' 

and the variance is given by 

(AL^)^Q / dt' dm' / dAm -^im' ,t'\m,t) — ; '-[L-iL)] . 

^ ' Jt-At Jo Jo dm'^ ' ' ' ^ dAmdt' \ ^ ') 

(11) 
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where Q is the normalisation factor (the compounded probabihty distribution in Eqs. (10) 
and (11) has been normahsed to 1. The effective lower limit for the integration over masses 
is set as described in Sect. (2.4) of CMT99. 

3. Results 

The results of our calculation are plotted in Fig. 1-3. 

In Fig. 1, we plot a direct comparison between the SSM, long-dashed line, with the 
MSSM, short dashed line, with the PEM ^, solid hne, and finally with the MPEM with a 
tilted CDM cosmogony. As well known, the SSM predicts that L ocT^ (Kaiser 1986), while 
the MSSM predicts non self-similar behaviour of the L-T relation: namely a L-T relation 
L (X ai scale of groups, L (x for rich clusters in agreement with observations and 
the L-T relation saturates toward L oc for higher temperatures. The plot shows that 
the MSSM predicts a similar behaviour of the L-T relation to that predicted by the PEM. 
Differences of maximum 10% are noted for smaller values of the temperature. 

It is noticed above, that a self-similar evolution for all clusters , at typical cluster 
temperatures (T > 2 keV), should lead to L ex T^, since free- free emission dominates the 
cooling. Instead the observed L-T relation is more steep, L oc 2"2.6-2.9^ meaning that lower 
temperature clusters and groups of galaxies are far less luminous than expected. Several 
different models have been proposed in order to explain the quoted behaviour in the L-T 
relation. The key-point of these models is that the X-ray luminosities of low-temperature 
clusters are small because their gas is less centrally concentrated than in hotter clusters, 
an effect that has been attributed to an universal minimum entropy level in intra-cluster 
gas, resulting from supernova heating (Ponman et al. 1999; Wu et al. 2000; CMT99), from 
heating by active nuclei (Wu et al. 2000) or from radiative cooling (Wu et al. 2000; Bryan 
2001; Pearce et al. 2000). In other terms for some reasons the core gas is less high than that 



^Thc PEM is based on hierarchical clustering. Group and cluster formation is envisaged in terms of 
DM potential wells evolving hierarchically, and engulfing outer baryons by accretion of smooth gas or by 
merging with other clumps. After a merging episode, the ICP in the wells falls back to a new, approximate 
hydrostatic equilibrium. This sequence of hydrostatic equilibria of the ICP is physically motivated for all 
merging events except for those involving comparable clumps (a mass ratio larger than ^ 1/4). However 
these sum up to less than 10% in the number. In the PEM, thermal energy of the infalling gas is initially 
due to stellar pre- heating (of nuclear origin); then it is increased to the virial value (of gravitational origin) 
when the accreted gas is bound in DM sub-clumps. So the prc-heating sets an effective threshold kTi ^ 0.5 
keV to gas inclusion, which breaks the self-similar correlation L ocT^ not only in its vicinity but also up to 
a few keV. 
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kT (keV) 



Fig. 1. — Comparison between the SSM, long-dashed hne, with the MSSM, short dashed 
hne, the PEM sohd hne, and MPEM dotted hne. 



kT (keV) 

Fig. 2. — MPEM model: the average L-T correlation with the 2 a dispersion (dotted lines), 
for a tilted cosmogony. Group data from Ponman et al. (1996) are represented by solid 
squares while cluster data from Markevitch (1998) are represented by sohd hexagons. 
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expected in the self-similar model. For example, an early episode of uniformly distributed 

supernova feedback, could rectify the problem by heating the un-condensed gas and therefore 
making it harder to compress in the core. In other words, the models with pre-heating and 
similar gives rise to the quoted break because they change the density in the core. During 
the hierarchical buildup an energy input pre heats the gas before it falls into new groups 
and clusters, so hindering its flow into the latter. The core density shall decrease and so the 
luminosity. 

A similar mechanism acts in the model of this paper. In fact, as shown in Del Popolo 
& Gambcra (1998), the angular momentum acquired by a shell centred on a peak in the 
CDM density distribution is anti-correlated with density: high-density peaks acquire less 
angular momentum than low-density peaks (Hoffman 1986; Ryden 1988). A greater amount 
of angular momentum acquired by low-density peaks (with respect to the high-density ones) 
implies that these peaks can more easily resist gravitational collapse and consequently it is 
more difficult for them to form structure, and in some conditions the structure formation by 
low mass peaks is even inhibited^^. 

The brecLk of the self-similcirity of the L-T relation may have also important 
consequences for determining the cluster masses from their luminosity. As shown 

by Shimizu et al. (2003), the predicted L-T relation is very sensitive to the assumed M-T 
relation, and then the non self-similarity of the L-T relation is strictly connected to that in 
the M-T relation. Also the M-T relation, as previously discussed, is non self-similar, and 
this behaviour has been interpreted in different ways (sec introduction). 
In Del Popolo 2002b, the bent in the M-T relation is entirely justified in terms of cluster 
tidal interaction with the neighbouring ones, or in other terms it is strictly connected to the 
asphericity of clusters (see Del Popolo & Gambera 1999 for a discussion on the relation be- 
tween angular momentum acquisition, asphericity and structure formation). Non-sphericity 
introduces an asymmetric bent, dependent on mass, in the M-T relation that gives rise to a 
different slope at the low mass end {T ~ 3keV): the lower the mass the larger the bent. 

The previous result is in agreement with AG result. In that paper, the authors used 
a nearly spherical collapsing region to obtain the M-T relation. According to their results, 
non-sphericity introduces an asymmetric, mass- dependent scatter (the lower the mass, the 
larger the scatter) for the M-T relation, thus altering the slope at the low masses end (T ~ 3 
keV). 



^''Onc interesting point to mention, at this point, is that several different assumptions are able to reproduce 
the observed L-T. This could mean that L-T relation, is not a very sensitive test: since almost any change 
to the "pure" self-similar model, reproduces this relation. 
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As commented in the introduction, heating/cooling mechanisms cire not nec- 
essciry to explain the observational L-T correlation. Bryan &: Norman (1998), 
carrying out large hydro dynamical simulations that followed the hierarchical 
evolution of clusters of galaxies, found that the observed M-T, L-T relations can 
be thoroughly reproduced if the number of particles and the spatial resolution 
are large enough. The importance of the numerical accuracy proves to be crucial 
to determining those relations, and so, for example, Mwanwong et al. (2001) 
using simulations with only 160^ peirticles and a spatial resolution of lOOh"^ Mpc, 
obtain that L ocT^ independent of mass, if radiative cooling is not implemented. 
In contrcist, Bryan &c Norman (1998), using simulations with 512^ particles and 
a spatial resolution of 50h~^ Mpc, did reproduce the observational bend of the 
L-T relation at the low-mass region. 

A priori, it is uncleeir why the L-T and M-T relations Eire that sensitive to 
resolution. One possible explanation goes in the direction of the results shown in 
this paper. Taking into account that proto-structures gain angular momentum 
owing to tidal interactions with other non-spherical structures, low resolution 
may hinder the gain of angular momentum by preventing an accurate deter- 
mination of the proto-structure shape. To clarify this point, we shall use the 
simulations of Mwanwong et al. (2001) and Bryan &c Norman (1998) as an ex- 
ample. The peirticle mass, rup, of the first was rup — 2.1 x lO^°/i~^M0, whereas the 
later used ~ 6 x 10^/i~^Mo in their simulations. Since a common chciracteristic 
of hydro-dynamical calculations is that all peirticles have the same maiss, clusters 
with kT < 2 keV (where the bend in the L-T relation starts to departure clecirly 
from self-similarity) have masses of M ~ 3 x lO^^hr^ Mq, i.e, they contain approx- 
imately < 1.4 X 10'' and < 5 x 10^ particles in Mwanwong et al. (2001) and 
Bryan Norman (1998) simulations, respectively. As we go to lower masses, we 
reduce the number of particles enclosed in proto-structures. As a consequence, 
the shape and, therefore, the inertia axes, may fluctuate randomly, which w 
ould lead in average to a systematical decrease of angulair momentum gained by 
low-mass proto-clusters. 

Besides the poorly determined shape of low-mass structures, one must also 
teike into account the possible effects of force resolution. In a typical N-Body 
evolution code, like for example the Treecode of L. Hernquist (1987), the force acting on a 
particle is given by the sum of two components: the force coming from the nearest neighbours 
and that coming from an expansion of the gravitational potential of the entire system up 
to quadrupole terms. As can be shown, the value of the average stochastic force in the 
simulation, Fsim, is an order of magnitude bigger than that obtained from the theory, Fth, of 
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stochastic forces. As a consequence only the higher force are taken into account, while the 
small fluctuations induced by the small scale substructure are not "seen". This is the case 
for CDM models in which the stochastic force generators are substructures at least three 
orders of magnitude smaller in size than the proto-structures in which they are embedded 
(e.g. clusters of galaxies). 

Taking into account the large simulations required to obtain a good descrip- 
tion of the L-T relation, it is not surprising that similar disagreements are reported in 
other cases too. For example, in the case of the M-T relation, it is noted that the results from 
different observational methods of mass measurements are not consistent with one another 
and with the simulation results (e.g., Horner, Mushotzky, & Scharf 1999, hereafter HMS; 
Neumann, & Arnaud 1999; Nevalainen, Markevitch, & Forman 2000, FRB). In general, X- 
ray mass estimates are about 80% lower than the predictions of hydro-simulations. On the 
other hand. X-ray mass estimates lead to normalisations about 50% higher than our result 
and simulations. 

One possible source for difference between theoretical and observational normalisations 
is that the values for are different in the two cases due to systematic selection effects. 
Also, intriguingly, Bryan & Norman (1998) showed that there is a systematic increase in the 
obtained value of (5 by increasing the resolution of the simulations. 

We would like to stress that even if the effects of angular momentum are not taken into 
account, this last process gives rise to self-similar structures only in a first approximation. 
In fact: (a) the effective spectral index rieff of CDM models depends, even if weakly, from the 
scale, going from values of rieff — —1.2 for clusters to ^efr — —2 for galaxies; (b) we live in a 
universe with cosmological constant different from zero, which means that there is a typical 
redshift at which it became important for cosmic dynamics; (c) even dark matter profiles 
are not perfectly self-similar, since they depends on the concentration parameter, which in 
turn is inversely proportional with mass, this because smaller structures formed, on average, 
at earlier times, when cosmic density was larger. 

As reported, Fig. 1 shows a slight difference between the SMMS prediction and that of 
PEM, being the slope predicted by SMSS, at low temperatures less steep than that of PEM, 
and MPEM. The difference is not so large, implying a difference in luminosity of 10% (larger 
for SMSS with respect to PEM). Also in Fig. 1, it is plotted the L-T relation predicted 
by MPEM. In this last case, the bending is produced by two effects: the threshold effect 
of the pre-heating temperature kTi ~ 0.5 keV (as in CMT99) and to the effect of angular 



^^/? = f5[l + f{l/P — l)f2b/f^m]) where / is the fraction of the baryonic matter in the hot gas, and fib is 
the density parameter of the baryonic matter. 
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momentum acquired by clusters. As a consequence, if we compare MPEM with PEM (or 
MSSM), the bending is larger (beside the threshold effect we have the acquisition of angular 
momentum) . 

Relative to this last item, looking at Fig. 1, one can see that the curve obtained from 
MSSM is very different to that corresponding to SSM, whereas the one from MPEM differs 
not much to that of PEM. The reason is the following: if we consider a cluster, without 
implementing pre-heating, the angular momentum acquisition is responsible for the slowing 
down and eventual stopping of the matter collapse towards the centre of the cluster, leading 
to the discussed consequences. Implementing prc-hcating, this gives rise (by heating the 
uncondensed gas and therefore making it harder to compress in the core) to a region at higher 
temperature and pressure that acts like a boundary for the infalling gas which, therefore, 
reduces the effects induced by angular momentum acquisition. 

In Fig. 2, we plot the results for the MPEM model: the average L-T correlation with the 
2 a dispersion (dotted lines), for a tilted cosmogony. Group data from Ponman et al. (1996) 
are represented by solid squares while cluster data from Markevitch (1998) are represented 
by solid hexagons. The L-T correlation is given by the double convolution (Eq. (10)), while 
AL is obtained by Eq. (11). The normahsation has been fitted on the data (see CMT99). 

The quantities and profiles of the PEM model are the same of CMT99, namely the 
reference cluster has a mass m — M/Mq, and a dark matter potential 0(r) as described in 
Appendix B. The density and temperature profiles are given by Eq. (30) and they should 
match the shock boundary conditions at the position r2 — itlyir. The average value and 
scatter of the parameter /3, given by Eq. (31), calculated through the PEM and shown in 
Fig. 2 of CMT99, increases from f3 = 0.5 to /5 ^ 0.9 while the baryonic fraction /2 is the 
one in Fig. 3 of CMT99. The 7 parameter is fixed as described in Appendix B. As the plot 
shows, in agreement with CMT99, the correlation is not a simple power law, but it starts 
as L oc for very rich clusters, and after it bends down with decreasing T. As previously 
told, the bending is induced by two mechanisms: the threshold imposed by the pre- 
heating temperature kTi ~ 0.5 keV (as in CMT99) and the angulcir momentum 
acquired by clusters. As a consequence, if we compare MSSM with PEM, the bending is 
larger (beside the threshold effect we have the acquisition of angular momentum). 

We want to point a similitude between the role of pre-heating temperature, T, in the 
PEM and that of the angular momentum £ in our model. In the PEM, the thermal energy 
of infaUing gas comes initially from stellar pre-heating (of nuclear origin) ; then it is increased 
to the virial value (of gravitational origin) when the accreted gas is bound to DM sub-clumps. 
So the pre-heating sets an effective threshold kTi ~ 0.5 keV to gas inclusion, which breaks 
the self-similar correlation L (xT^ not only in its vicinity but also up to a few keV. Increasing 
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the pre-heating temperature the bending in the L-T relation becomes more pronounced. A 
similar process occurs if the acquired angular momentum is larger. 

A similar fitting formula to that of CMT98, for the predicted L-T correlation (for 
T > Ti ~ 1 keV), is given by: 

L = a^T^+^Mp/Po)'/' (12) 
aL oc fi°-3(l + ^)0-22M' + (l-Qo)e-°-'(^+^) 
aL = ai(l + z)-°-^e-«^(^-^i)/"S-^(^+^)"\ 

where the luminosity is expressed in units of 10^^ erg/s and the temperature in keV, and 
with oi = 1.2 and 02 = 0.17. 

At temperatures larger than the threshold kTi ~ 0.5 keV, the relative AL/L remains 
constant around 25%. A study of the dependence of (L) and AL on shows that both 
these quantities increase with increasing VLq, similarly to what shown in CMT99 

Fig. 3, shows MSSM compared with observational data. Similarly to Fig. 2, we plot 
the average L-T correlation with the 2 a dispersion (dotted lines), for a tilted cosmogony. 
Group data from Ponman et al. (1996) are represented by sohd squares while cluster data 
from Markevitch (1998) are represented by sohd hexagons. The L-T correlation can be fitted 
in this case by a similar formula to that of Eq. (13), with oi = 1.28 and 02 = 0.19. As 
reported. Fig. 1 shows a slight difference between the SMMS prediction and that of PEM, 
being the slope predicted by SMSS, at low temperatures less steep than that of PEM. The 
difference is not so large, implying a difference in luminosity of 10% (larger for SMSS). The 
fit to data of SMSS model, as Fig. 3 shows, is also very good. 

To summcirise, the key idea of the SMSS model and of other mechanisms 
proposed to reproduce the non-self-similarity of the L-T relation is in all Ccises 
fairly similar: If one wants to have clusters less luminous than SSM prediction, it is neces- 
sary to have a physical process that reduces the quantity of gas infalling towards the centre of 
the cluster which, therefore, reduces the core luminosity. In the case of heating/cooling 
models, some energy input pre-heats the gas before it falls into new groups and 
clusters, hindering its flow into the latter. In the SMSS model, that role is played 
by the initial spin present in proto-clusters. 



^^This is because the underlying strength of the current shocks grows on average as the merging rate 
(moderately) increases on approaching the critical cosmology, see Lacey & Cole (1993). 
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4. Conclusions 

In this paper we showed that the presence of angular momentum during the collapse of 
a proto-structure leads to a non-self-similar L-T relation. The quoted effect leads, in X-rays, 
to a luminosity-temperature relation that scales as L oc T^, at scale of groups, flattening to 
L ocT^ for rich clusters and converging to L oc at higher temperatures. 

These results are in disagreement with the largely accepted assumption that heat- 
ing/cooling processes and similar are fundamental in the originating the non-self similar 
behaviour (shaping) of the L-T relation. As Bryan &: Norman (1998) showed, it is not 
necesseiry to hypothesise pre-heating/cooling models in order to reproduce ob- 
servations, on the contrciry, it is possible to reproduce the observed L-T relation 
if the spatial and mass resolution Eire accurate enough. Poorly resolved clusters, 
with few pEirticles enclosed, lead to self-similcir L-T curves. 

We have shown that the large bend of the L-T relation is caused by the fact that the 
angular momentum acquired by a shell centred on a peak in the CDM density distribution is 
anti-correlated with density: high-density peaks acquire less angular momentum than low- 
density peaks. A greater amount of angular momentum acquired by low-density peaks (with 
respect to the high-density ones) implies that these peaks can more easily resist gravitational 
collapse and consequently it is more difficult for them to form structure. This results in a 
tendency for less dense regions to accrete less mass with respect to a classical spherical model. 
As a consequence, the X-ray luminosities of low-temperature clusters are small because their 
gas is less centrally concentrated than in hotter clusters. 
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5. Appendix A: M-T relation 

As previously quoted, numerical methods and simple scaling arguments suggest that 
the X-ray temperature of clusters, Tx, can be directly related to their masses as Mvir oc 
Tx^^Pi^^^^A~jy^, where pb is the critical density, A^ir the mean density within the virial 
radius i?vir- 

In Del Popolo & Gambera (1999) and Del Popolo (2002), we got the M-T relation in 
two different ways: (1) modifying the top-hat model; (2) modifying Voit & Donahue (1998) 
(hereafter V98) model. 

In the first, we modified the top-hat model in order to take account of angular momentum 
acquisition by protostructures and used a modified version of the virial theorem in order 
to include a surface pressure term (V2000, AC). This correction is due to the fact that at 
the virial radius i?vir the density is non-zero and this requires a surface pressure term to be 
included in the virial theorem (Carlberg, Yee & Ellingson 1997) (the existence of this con- 
fining pressure is usually not accounted for in the top- hat collapse model). The derivation of 
the previous relation is fundamentally based on the approximation of cluster formation with 
the evolution of a spherical top-hat density perturbation (Peebles 1993) and on the addi- 
tional assumption that each cluster observed at a redshift z has just reached the moment of 
virialization. This last assumption is currently known as the late-formation approximation, 
which is a good one in a critical = because for this value of massive clusters develop 
rapidly at all redshifts and the moment of virialization is always close to that of observation. 
In other terms for f2o = 1, the accretion rate remains sufficiently high, and this implies that 
the clusters we actually observe attained their observed masses recently. In the i^o < 1 case 
cluster formation is "shutting down" and it is necessary to take account of the differences 
between the moment of virialization and that of observation. The problem becomes worse 
going through f2o << 1: in fact in the late- formation approximation Myir rises steadily since 
PbAvir declines indefinitely, while we expect that the cluster formation is going to stop 

The late-formation approximation is a good one for many purposes, but a better one 
can be obtained in the low-Q limit. As can be found in the literature, there are two ways 
of improving the quoted model. One is to define a formation redshift Zi at which a cluster 
virializes and after the properties of observed clusters at z are obtained by integrating over 
the appropriate distribution of formation redshifts (Kitayama & Suto 1996; Viana h Liddle 
1996). The second possibility is the one described by V98, V2000. In this approach, the top- 
hat cluster formation model is substituted by a model of cluster formation from spherically 



The result of the late-formation approximation is displayed in Eqs. 18-19 of Del Popolo (2002) 
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symmetric perturbations with negative radial density gradients. The fact that clusters form 
gradually, and not instantaneously, is taken into account in the merging-halo formalism of 
Lacey & Cole (1993). In hierarchical models for structure formation, the growth of the largest 
clusters is quasi-continuous since these large objects are so rare that they almost never merge 
with another cluster of similar size (Lacey & Cole 1993). So, Lacey & Cole (1993) approach 
extends the Press-Schechter formalism by considering how clusters grow via accretion of 
smaller virialized objects. Summarising, in order to obtain the proper normalisation and 
time evolution of the M-T relation, one has to account: 

a) for the continuous accretion of mass of clusters; 

b) for the non-zero density at i?vir, requiring a change in the virial theorem by including a 
surface pressure term. 

The M-T relation derived by means of a model of continuous accretion, differs from the 
late-formation model in both normalisation and time-dependent behaviour. 

In order to obtain the M-T relation in this second approach, we assume, as shown by 
V2000, the mass grows hke M oc u;-3/("+3) (Lacey & Cole 1993; V98; V2000). The virial 
energy of the cluster, —E, can be calculated by integrating the 

In order to obtain an expression for the kinetic energy, we firstly calculated E/M: 

(13) 



M M 10(mi - 1) V 



E f edM 3mi /27rG\ - ,^2 



mi 



2 



t y M8/3 



Ki{mi,x) = (mi — 1) FxLerchPhi{x, 1, 3mi/5 + 1) — 

(mi - 1) FLerchPhi{x, 1, 3mi/5) (14) 



^^A comparison of the normalisation predicted by the late-formation model with that predicted by simu- 
lations of Evrard, Mctzlcr & Navarro (1996) shows that when Oq = 1 this normalisation is only 4% below 
the empirical value, but it lies 20% below it for fig = 0.2. In the case of V2000 model and for a power-law 
spectrum, a comparison with the same simulations show that the temperature normalisation of the n = —2 
case deviates by less than 10% over the range 0.2 < Oq < 1 and by ~ 18% in the case n = — 1 (V2000). 
The normalisation obtained by the V2000 model, even if it is more accurate than that given by the late for- 
mation, or that by AC which is in agreement with hydro-simulations, show a noteworthy discrepancy when 
compared with X-ray mass estimates (about 50% for the AC model; see also V2000). One possible source 
for differences in theoretical and observational normalisations may be due to the fact that /3 is different in 
the two cases because of systematic selection effects. For example, as shown by Bryan & Norman (1998), 
increasing the resolution of simulations there is an increase in the value of /3. So summarising, for what 
concerns normalisation, the continuous formation model gives more precise results than the late formation 
one, but in any case if we want to fit observations we need to shift the normalisation (see AC). 
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where, 

nil = 5/(n + 3), tn = ^ 3-, x = 1 + which is connected to mass by M = 

MoX~^"*i/^ (V2000) and ^ — where rta is the turn-around radius and xi is defined by the 

relation M = with P6 is the background density. The LerchPhi function is defined as 

follows: 

LerchPhi{z, a,v)^Yl u + ^lY ^^^^ 

n=0 ^ ' 

The angular momentum, >C, acquired by the protostructure, is calculated using the same 
model (and same spectrum) as described in Del Popolo & Gambera (1998, 1999). More 
hints on the model and some of the model limits can be found in Del Popolo, Ercan & 
Gambera (2001)). Then ,the virial theorem with the surface pressure term correction, as in 
V2000, is used in order to get a connection between the kinetic energy and temperature. We 
utilise the usual relation: 

(K) = ^ (17) 

(AC), where k is the Boltzmann constant, /i — 0.59 is the mean molecular weight, rUp the 
proton mass and f3 — ^^J^^ , being ay the mass- weighted mean velocity dispersion of dark 
matter particles, and P = /?[1 + /(1//5 — l)fli,/flm], where / is the fraction of the baryonic 
matter in the hot gas, and f^b is the density parameter of the baryonic matter. In this way, 
we finally get: 



5 2P mi-l\ tn J 



1 ^ f^Jl]^ + Ki{mi,x) 



mi V ^ / Mq^^ 



(18) 



where a = 2p{R- )~p ratio between kinetic and total energy (V2000). If Ki — 0, Eq. 

(13) reduces to Eq. (10) of V2000. As stressed by V2000, some factors give rise to an higher 
value of E/M with respect the case of the late-formation value. The mi/(mi — 1) value 
which accounts for the effect of early infall. The 1/mi value in the square bracket of Eq. 
(13) which accounts for the cessation of cluster formation when t » to,. Finally in Eq. (13) 
a new term is present, which comes from the tidal interaction. 

Using the relation Ayir = (see V2000), and in the early-time limit: {t « ta), Eq. 
(18), reduces to: 



which, in the case n ~ — 2, a ~ 2 is identical to the late- formation formula, described in 
V2000 (see their Eq. (8)). Normalising Eq. (18) similarly to V2000, we get: 



kT ~ ^keV 
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(20) 



where KQ{mi,x) indicates that Ki{mi, x) must be calculated assuming t — to 

Eq. (20) when compared to the result of V2000 (Eq. 17) shows an additional term, 
mass dependent. This means that, as in the case of the top-hat model, the M-T relation is 
no longer self-similar showing a break at the low mass end (see next section). 



The X-ray bolometric luminosity of a cluster is given by Eq. (1), which in CMT98 
notation is: 



Here T(r) is temperature in the plasma and r2 is the cluster boundary, that we take to 
be close to the virial radius -Rvir oc M^/^p^^/'^, where p{z) oc (1 + zY is the DM density 
in the cluster, proportional to the average cosmic DM density Pu{z) at formation. The 
infalling gas is expected to become supersonic near r2 (see, e.g., Perrenod 1980; Takizawa & 
Mineshige 1998) so that a shock front will form there. The conservations across the shock of 
mass, energy and stresses yield the Rankine-Hugoniot conditions, i.e., the temperature and 
density jumps from the outer values Ti and rii to T2 and ^2 just interior to r2- Then the 
luminosity may be rewritten in the form 



where x = r/r2. rii is fixed by rii oc fuPu/f^p, in terms of the universal baryonic fraction 
fu] whereas Ti is determined only statistically, through the diverse merging histories ending 
up in the mass M. In sum, a given dark mass M admits a set of ICP equilibrium states 
characterized by different boundary conditions, each corresponding to a different realization 
of the dynamical merging history. It is the convolution over such set which provides the 
average values of L and Rx, and their scatter. Following CMT98, the pre-shock temperature 
in a merging event is that of the infalling gas, and if the latter is contained in a sufficient deep 
potential well, Ti is the virial temperature Tiy oc Am/r of the secondary merging partner; 
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(21) 




(22) 
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on using r oc {Am/pY^^ this writes 

kTi^ = 4.5 (Am)2/3 (p/p„)V3 keV, (23) 

where the numerical coefficient is taken from Hjorth, Oukbir & van Kampen (1998), and the 
masses m — M/Mq are normahsed to the current value Mq = 0.6 x IO^^Qq^"^-^© (i-e., the 
mass enclosed within a sphere of 8/i~^ Mpc), so in the following the actual value of Ti will 
be 

Ti = maa;[Ti^,ri,]i5 (24) 

Given Ti, the boundary conditions for the ICP in the cluster is set by the strength of the 
shocks separating the inner from the infalling gas. In the case of three degrees of freedom 
and for a nearly hydrostatic post-shock condition with V2 «Vi, assuming the shock velocity 
to match the growth rate of the virial radius R^i^{t): 

r(l + v^lTi)^ 7 3 6^ 1 

- y 1 + 1^' " ^ (TTTTTlj^J ■ ^ ^ 

Cavaliere, Menci e Tozzi (1997) (CMT97). 

Here e = 15kTi/A/jmpvf and /j is the average molecular weight; the inflow velocity Vi 
is set by the potential drop across the region of nearly free fall, to read Vi ~ a/ —<t>2/fnp in 
terms of the potential 4>2 at r^- In the case of strong shocks, appropriate to "cold inflow", 
e << 1, as in rich clusters accreting small clumps and diffuse gas, the approximation 

kT2^-(j)2/3 + 3kTi/2 (26) 

holds, where 02 is the gravitational potential energy at r2 — -Rvir- For e > 1 the shock is 
weak, and T2 — T^. Prom T2 and Ti, the density jump at the boundary n2/ni is found to 
read (see CMT97) 

!!a.2(l-^)+^4(l-^)^^l^^ (27) 

m V T2) V \ T2J TsJ ^ ' 

Adopting the polytropic temperature description T{x)/T2 = [n{x) / n2\^~^ , with the 
index 7 in the range 1 < 7 < 5/3, and that the radius r2 may be written in terms of 
temperature oc m/r2 and that m oc pr|, leading to r2 oc {t/pY^'^, the luminosity can be 
written in the form: 



^'"^An independent lower bound kTi^ w 0.5 keV is provided by preheating of diffuse external gas, due to 
feedback energy inputs following star formation and evolution all the way to supernovae (David et al. 1995; 
Renzini 1997). 
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L oc 



n2 

ni 



m 



P 



To 



[n{r)/n2] 



2+(7-l)/2 



(28) 



where the bar denotes the integration over the emitting volume < and p is the average 
DM density in the cluster, proportional to and so to rii. 

Eq. (28) can be also cast in the form: 



L OC ( — p 



T 



1/2. 



[n{r)/n2 



|2+(7-l)/2^4/3 



K 



S73 



(29) 



The ratio n{x)/n2 is obtained from the hydrostatic equilibrium dP/rripndr = —GM{< 
r)/r'^ = ^d(t)/dr with the polytropic pressure P(r) = kT2n2[n{r) /n2\' . This yields (see 
Cavaliere & Fusco Femiano 1978; Sarazin 1988, and bibliography therein) the profiles 



n{r) 
n2 



\T(r)^ 



1/(7-1) 



{l + ^/5[02-0(r)]} 



1/(7-1) 



(30) 



where (j) = (p/pmpal is the potential normalised to the associated one-dimensional DM 
velocity dispersion at r2- The ICP disposition in eq. (11) relative to the DM depends on the 
parameter, already met previously: 



(3 = pmpa2/ kT2 , 



(31) 



and is further modulated by the second parameter 7, to yield as the latter increases flatter 
profiles n(r) and steeper T(r). 

The function /3(T) can be easily computed from Eq. (26) for a given dark matter 
potential 02 corresponding to p(r). 0(r) and a(r) are obtained in agreement with Navarro, 
Prenk & White (1997). 



^^For the King potential (see Sarazin 1988) and CMT97, with core radius rc = i?„/12, /3(T) increases from 
/3 ~ 0.5, for T ~ Ti to /3 ~ 0.9 for T » Ti. A similar result is obtained for a Navarro et al (1996) potential. 
The other parameter 7 will be bounded according to CMT99, The polytropic index 7 > 1 describes the 
equation of state for the ICP. An upper bound to it arises if the overall thermal energy of the ICP is not to 
exceed its gravitational energy. The thermal and the gravitational energy are computed using the profiles in 
Eq. (30), and their ratio is given in Fig. 4 of CMT99, to show that the upper bound 7 < 1.3 holds. It turns 
out that observations by Markevitch et al. (1997) are consistent with the T(r) predicted when 7 = 1.2 ±0.1, 
in our allowed range. Hereafter we shall focus on 7 = 1.2. 
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7. Appendix C: Calculation of the anguleir momentum 

The effect of tidal torques on structures evolution has been studied in several papers 
especially in connection with the origin of galaxies rotation (Hoyle 1949; Peebles 1969; White 
1984; Ryden 1988 (hereafter R88); Eisenstein & Loeb 1995). 

Following Eisenstein & Loeb (1995), we separate the universe into two disjoint parts: 
the collapsing region, characterised by having high density, and the rest of the universe. The 
boundary between these two regions is taken to be a sphere centred on the origin. As usual, 
in the following, we denote with p(x), being x the position vector, the density as function 
of space and 5(x) = The gravitational force exerted on the spherical central region 

by the external universe can be calculated by expanding the potential, $(x), in spherical 
harmonics. Assuming that the sphere has radius R, we have: 

oo , / 

^(^) = E ^ E ^M^UO: <t>)x' (32) 

l=Q m=-l 

where Yi^ are spherical harmonics and the tidal moments, aim, are given by: 

poo 

aUx)^p^ / YUe,(l>)p{s)s-'-'dh (33) 
Jr 

In this approach the proto-structure is divided into a series of mass shells and the torque 
on each mass shell is computed separately. The density profile of each proto-structure is ap- 
proximated by the superposition of a spherical profile, S{r). and a random CDM distribution, 
e(r), which provides the quadrupole moment of the proto-structure. To the first order, the 
initial density can be represented by: 

p(r)=pb[l + 5(r)][l+£(r)] (34) 

where £(r) is given by: 

(k^r) = m (35) 

being P{k) the power spectrum. The torque on a thin spherical shell of internal radius x is 
given by: 

t{x) = ^ / £(x)xxv$(x)cin (36) 

where Mgh — 47rpb [1 + ^{x)] x^Sx. Before going on, I want to recall that we are interested 
in the acquisition of angular momentum from the inner region, and for this purpose we take 
account only of the I — 2 (quadrupole) term. In fact, the I — term produces no force, while 
the dipole {I — 1) cannot change the shape or induce any rotation of the inner region. As 
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shown by Eisenstein & Loeb (1995), in the standard CDM scenario the dipole is generated 
at large scales, so the object we are studying and its neighbourhood move as bulk flow with 
the consequence that the angular distribution of matter will be very small, then the dipole 
terms can be ignored. Because of the isotropy of the random field, s{x), Equation (36) can 
be written as: 

< kr >= a/(30)^ [< a2m{xf >< q2m{xf >-< a2m{x)qlm{x) >'^f'^ (37) 

where <> indicates a mean value of the physical quantity considered. As stressed in the 
next section, following Eisenstein & Loeb (1995), the integration of the equations of motion 
shall be ended at some time before the inner external tidal shell (i.e., the innermost shell of 
the part of the universe outside the sphere containing the ellipsoid) collapses. Then the inner 
region behaves as a density peak. This last point is an important one in the development of 
the present paper. 

An important question to ask, before going on, regards the role of triaxiality of the el- 
lipsoid (density peak) in generating a quadrupole moment. Equation (37) takes into account 
the quadrupole moment coming from the secondary perturbation near the peak. The den- 
sity distribution around the inner region is characterised by a mean spherical distribution, 
5, and a random isotropic field. In reality the central region is a triaxial ellipsoid. It is then 
important to evaluate the contribution to the quadrupole moment due to the triaxiality. 
Remembering that the quadrupole moments are given by: 

?2n. = / Y^e, 4>)s'p{s)dh = ^ / y;^(^, m^)dn (38) 
J|r|<R 47r J 

and approximating the density profile as: 

5(x) =< 5{x) > Spherical f {x)A{e,p) (39) 

being < 5{x) > spherical thc mean spherical profile, i^ = ^ the peak height and a the r.m.s. 
value of 5. The function A[e,p) of the triaxiality parameters, e and p, is given by: 

A{e, p) = 3e(l - sin^ 9 - sin^ sin^ 0) + p(l - 3 sin^ cos^ (f)) (40) 

while the function f{x) is given (R88) by: 

where ^, a and i?* are respectively the two-point correlation function, the mass variance and 
a parameter connected to the spectral moments (see Bardeen et al. 1986, equation (4.6d), 
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hereafter BBKS). Substituting equation (39) and equation (40) in equation (38) it is easy to 
show that the sum of the mean quadrupole moments due to triaxiahty is: 

^ ^ < q^^ix) >= vx^f{x) (^i-v^(e -p) + l-^/^{Ze + p)^ (42) 

which must be compared with that produced by the secondary perturbations, e: 

< q2m{xr J k'P{k)n{kxfdk (43) 

where j2 is the Bessel function of order 2. The values of e and p can be obtained from the 
distribution of elhpticity and prolateness (BBKS, equation (7.6) and figure 7) or for i/ > 2 
by: 

e = ^- T7^ (44) 

v^x [1 + 6/(5x2)]^/' 

and 

p 2 (45) 

5x4 [1 + 6/(5x2)]2 V ) 

(BBKS equation (7.7)). where x is given in BBKS (equation (6.13)). In the case of a peak 
with = 3, wc have e ^ 0.15, p ~ 0.014 while for peaks having u = 2 and v = 1 they are 
respectively given by e ~ 0.2, p ~ 0.03 and e ~ 0.25 p ~ 0.04. 

As shown in figure 1 of Del Popolo et al. (2001), for a 3(7 profile, the source of quadrupole 
moment due to triaxiahty is less important than that produced by the random perturbations 
e in all the proto-structure, except in the central regions where the quadrupole moment due 
to triaxiahty is comparable in magnitude to that due to secondary perturbations. In other 
words, the triaxiahty has a significant effect only in the very central regions, which contains 
no more than a few percent of the total mass and where the acquisition of angular momentum 
is negligible. It follows that the triaxiahty can be ignored while computing both expansion 
and spin growth (R88). Moreover, as observed by Eisentein & Loeb (1995), the ellipsoid 
model docs better in describing low shear regions (having higher values of z/), whose collapse 
is more spherical and then the effects of triaxiahty are less evident. Just this peaks, having 
at least z/ > 2, shall be studied in this paper. In any case, even if the triaxiahty was not 
negligible it should contribute to increment the acquisition of angular momentum (Eisenstein 
& Loeb 1995), and finally to a larger effect on the density evolution, (i.e., a larger reduction 
of the growing rate of the density) . 

In order to find the total angular momentum imparted to a mass shell by tidal torques, it 
is necessary to know the time dependence of the torque. This can be done connecting g2m and 
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a2m to parameters of the spherical collapse model (Eisenstein & Loeb 1995 (equation (32), 
R88 (equation (32) and (34)). Following R88 we have: 

/i(^) - (I) fM 



and 



4X4/3 



a2ra{0) = a2m,0 ( 3 ) 5o(^ " Siu e) ^ (47) 



The collapse parameter 9 is given by: 

t{e) = -toSo''^\e-sme) (48) 

Equation (46) and (47), by means of equation (37), give to us the tidal torque: 

where /i(^) and /2(^) are given in R88 (Eq. 31), tq and 5o = are respectively the torque 
and the mean fractional density excess inside the shell, as measured at current epoch to- The 
angular momentum acquired during expansion can then be obtained integrating the torque 
over time: 

^ = /t(^)|^^ (50) 

As remarked in the Del Popolo et al. (2001) the angular momentum obtained from equation(50) 
is evaluated at the time of maximum expansion tM- Then the calculation of the angular mo- 
mentum can be solved by means of equation (50), once we have made a choose for the power 
spectrum. With the power spectrum and the parameters given in the next section and for 
a = 2 peak, the model gives a value of 2.5 x lO^^gcm^/s. As previously quoted, we as- 
sume that from tu on, the ellipsoid has this constant angular momentum. Following the 
procedures 1) and/or 2), we shall be able to get the time evolution of the density. 



